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Abstract. We study the Boltzmann equation for a space- homogeneous gas of 
inelastic hard spheres, with a diffusive term representing a random background 
forcing. Under the assumption that the initial datum is a nonnegative L^(]R^) 
function, with bounded mass and kinetic energy (second moment), we prove the 
existence of a solution to this model, which instantaneously becomes smooth and 
rapidly decaying. Under a weak additional assumption of bounded third moment, 
the solution is shown to be unique. We also establish the existence (but not 
uniqueness) of a stationary solution. In addition we show that the high-velocity 
tails of both the stationary and time-dependent particle distribution functions 
are overpopulated with respect to the Maxwellian distribution, as conjectured by 
previous authors, and we prove pointwise lower estimates for the solutions. 



Introduction 

In recent years a significant interest has been focused on tlie study of kinetic mod- 
els for granular flows pUl ^] ■ Depending on the external conditions (geometry, 
gravity, interactions with surface of a vessel) granular systems may be in a variety 
of regimes, displaying typical features of solids, liquids or gases and also producing 
quite surprising effects |SHj. Finding a systematic way to describe such systems 
under different conditions is a physical problem of considerable importance. At the 
same time, recent developments in this area gave rise to several novel mathematical 
models with interesting properties. 

In the case of rapid, dilute flows, the binary collisions between particles may be 
considered the main mechanism of inter-particle interactions in the system. In such 
cases methods of the kinetic theory of rarefied gases, based on the Boltzmann- Enskog 
equations have been applied j211 EHl EOl • 

A very important feature of inter-particle interactions in granular flows is their in- 
elastic character: the total kinetic energy is generally not preserved in the collisions. 
Therefore, in order to keep the system out of the "freezing" state, when particles 
cease to move and the system becomes static, a certain driving mechanism, supply- 
ing the system with energy, is required. Physically realistic driven regimes include 
excitation from the moving boundary, through-flow of air, fluidized beds, gravity, 
and other special conditions. We accept a simple model for a driving mechanism. 
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the so-called thermal bath, in which we assume that the particles are subject to un- 
correlated random accelerations between the coUisions. Such a model was studied 
in 1301 one- dimensional case, and in [SZj in general dimension. 

We study the model [SZj in the space-homogeneous regime, described by the 
following equation: 

(0.1) dJ-fiAJ = Q{f,f), veR'', t>0. 

Here / is the one-particle distribution function (particle density function in the 
phase space), which is a nonnegative function of the microscopic velocity v and the 
time t] we shall assume N > 2 (dimension 1 could be treated as well but would 
require a few notational changes). On the right-hand side of equation (jO.lll there is 
the inelastic Boltzmann-Enskog operator for hard spheres (the details of which are 
given below); the term —fiA^f, /i = const, represents the effect of the heat bath. 
Without loss of generality we can set yU = 1 (see Section Fl.Sj) . which we will from 
now on assume. In the sequel, we shall often abbreviate into just A. 

One of the interesting features of the model (jO.lj) is the fact that it possesses 
nontrivial steady states described by the balance between the collisions and the 
thermal bath forcing. Such steady states are given by solutions of the equation 

(0.2) /iA,/ + Q(/,/) = o, vem"". 

Solutions of ()0.2j) have been studied in |37j by means of formal expansions. The 
same problem was also studied in [21 and in jH], for a different kind of interactions, 
namely the Maxwell pseudo-particle model 123 1211; niethods of expansions 
and the Fourier transforms, respectively. In reference fT] the rigorous existence of 
radially symmetric steady solutions for the Maxwell model was established. 

The aim of this study is to develop a rigorous theory of for the inelastic hard 
sphere model, and to investigate the regularity and qualitative properties of the 
solutions. We prove that equation has a unique weak solution under basic 

assumptions that the initial data have bounded mass and kinetic energy, and satisfy 
some additional conditions (bounded entropy for existence, L^(R^) for regularity, 
and bounded third moment in for uniqueness). The thermal bath (diffusion) term 
in ()0.1|l is responsible for the parabolic regularity of solutions: the weak solutions 
become smooth, classical solutions after arbitrarily short time. We apply generally 
similar techniques, based on elliptic regularity, to treat the steady case. Finally, we 
establish lower bounds, for both steady and time- dependent solutions, proving that 
the distribution tails are "overpopulated" with respect to the Maxwellian, as was 
suggested in [SI]. The lower bound for steady solutions is given by a "stretched 
exponential" y4exp(— a|f p/^), with a = a{a,ix). In the time-dependent case the 
bound holds with A = A{t), where A{t) is a generally decaying function of time. 
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We emphasize that the appearance of the "3/2" exponent is a specific feature 
of the hard sphere model with diffusion, and could be predicted by dimensional 
arguments (cf. [SZI)- On the other hand, the Maxwell model with diffusion results 
in a high-velocity tail with asymptotic behavior Cexp{—c\v\), see As a general 
rule, the exponents in the tails are expected to depend on the driving and collision 
mechanisms j2 CHI CZl El • In fact, deviations of the steady states of granular systems 
from Maxwellian equilibria ( "thickening of tails" ) is one of the characteristic features 
of dynamics of granular systems, and has been an object of intensive study in the 
recent years [23 EH E3 E21 ■ 

We remark that the "3/2" bound has rather important practical implications 
as well. In particular, it indicates that the approximate solutions based on the 
truncated expansion of the deviation from the Maxwellian into Sonine polynomials 
jSZl El IS2] could only be valid for moderate values of |fp. Any conclusions about 
the tail behavior drawn from such an expansion should be questioned. Indeed, since 
the deviation function is growing rapidly for |f | large (it is in the weighted space, 
but not in !), the Sonine polynomial expansion should in general be expected to 
have poor approximation properties in this region. 

The paper is organized as follows. The first section contains the preliminaries, 
where we introduce the inelastic collision operator and establish several basic iden- 
tities which are important in the sequel. In section 2 we establish the bounds for 
the energy and entropy of solutions In Section 3 we study the moments of the distri- 
bution function by analyzing the moment inequalities for equations ()0.1|) and ()0.2|) . 
The key point in analyzing the moments is the so-called Povzner inequalities, well- 
known for the classical Boltzmann equation jHH 1^ El IBHl IS IHUI , which we here 
extend to the case of inelastic interactions and present in a general setting of poly- 
nomially increasing convex test functions. In Section 4 we study the estimates of 
the inelastic collision operator in spaces with polynomial weights, extending the 
results in j2l] to the inelastic hard sphere case. We continue by establishing apriori 
regularity estimates, based on the interpolation of spaces and the Sobolev-type 
inequalities. In Section 5 we present a rigorous proof of the existence and regularity 
of the time-dependent and steady solutions. The arguments presented there also 
justify the formal manipulations performed in Sections 2, 3 and 4. In Section 6 we 
show the uniqueness for the time-dependent problem using Gronwall's lemma. Fi- 
nally, in Section 7 we compute lower bounds for the stationary and time-dependent 
solutions. 



1. Preliminaries 



1.1. Binary inelastic collisions. We study the dynamics of inelastic identical hard 
balls with the following law of interactions. Let v and be the velocities of two 
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particles before a collision, and denote by m = f — f * their relative velocity. Let the 
prime symbol denote the same quantities after the collision. Then we assume 

(u' ■ n) = —a (u ■ n), 

(1-1) / / 

u — {u ■ n) = u — [u- n), 

where n is the unit vector in the direction of impact, and < a < 1 is a con- 
stant called the coefficient of normal restitution. Setting w = v + v^, and using the 
momentum conservation we can express v' and v'^ as follows: 

, w u' , w u' 

(1.2) -y'=- + — , V'= , 

\ ) 2 2 * 2 2 

By substituting into ()1.2|) and equations the post-collisional velocities v' 

and v'^^ are uniquely determined by the pre-collisional ones, v and f*, and the impact 
parameter n (cf. jTUI, |SZ|)- 

The geometry of the inelastic collisions defined by relations (jl.lj) . ()1.2|) is shown in 
Figure 1. For every v and f * fixed, the sets of possible outcomes for post-collisional 
velocities are two (distinct) spheres of diameter Thus, it is convenient to 

parametrize the relative velocity after collision as follows: 

(1.3) u' = {I- (3)u + (3\u\a, 

where we denoted /3 = ^-t^. The relations p.2|) and ()1.3|) define the post-collisional 
velocities in terms of v, v^, and the angular parameter a G S^~^. 

1.2. Weak form of the collision operator. We define the collision operator by 
its action on test functions, or observables. Taking ip = 'ip{v,t) to be a suitably 
regular test function, we introduce the following weak bilinear form of the collision 
term: 



(1.4) / Q{g, f)ipdv = / / / fg^,{ip'-ip)\u\b{u,a)dadvdv^,. 

Here and below we use the shorthand notations / = f{v,t), = g{v^,,t), ip' = 
ip{v',t), etc. The function b{u,a) in ()1.6|1 is the product of the Enskog correlation 
factor k{p, d) (which is a constant in the space-homogeneous case) by the differential 
collision cross-section, expressed in the variables u, a. In the case of hard-sphere 
interactions, 



b{u, a) = k{p, d) I - 



2 J V 2 

where v = u/\u\, and d is the diameter of the particles. Notice that the hard 
sphere cross-section depends only on the angle between u and a, and is generally 
anisotropic, unless = 3. Without restricting generality, by choosing the value of 
d accordingly, we can always assume that 

(1.5) / b{u,cx)da=l. 
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All lengths shown in assumption \u\/2 = 1. Unit vectors n, a, uo not 
to scale. 



Of course, to write down the Boltzmann operator we only need /), but later 
on it will be sometimes convenient to work with the bilinear form Q{g,f). An 
explicit form of Q will be given later on; however for many purposes it will be easier 
to work with the weak formulation which is also quite natural from the physical 
point of view (it is analogous to the well-known Maxwell form of the Boltzmann 
collision operator jSHl Chapter 1, Section 2.3]). 

In the case when f = g in ()1.4|) . we can further symmetrize and write 

(1.6)/ Q{f,f)tpdv = l [ [ [ fM^'+^:-^-^,)\u\b{u,a)dadvdv,. 

Notice that the particular form of the inelastic collision laws enters ()l.fi|l only through 
the test function ip'. 

1.3. Equations for observables and conservation relations. Using the weak 
form ()1.6|) allows us to study equations for average values of observables given by 
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the functionals of the form J^pf ftp dv. Namely, multiplying equation ()0.1|) by a test 
function ipiv^t) and integrating by parts we obtain 



t=T 



T r rT 



1-7) / fijdv - / / f{dtilJ + ^.ij)dvdt= / / Q{f,f)i:dvdt. 



*=0 Jo JR^ Jo 

With the weak form of the collision operator, it is easy to verify formally 

the basic conservation relations that follow from ()0.1|) . Namely, setting ip = 1 and 
ip = Vi in ()1.7j) and assuming that f^^f f ip dv is differentiable in t, we obtain the 

conservation of mass and momentum: 

(1.8) ^ / f{l,vi,...,VN}dv = 0. 

at j^N 

Further, taking ip = |f p and computing 

n n^ I /|2 , I M2 I |2 I |2 1 - 1 - ('^ ■ 0") I ,2 
(1-9) \v\ + \vj -\v\ -\v*\ = |m| , 

we obtain the following relation for the dissipation of kinetic energy: 
(1.10) ^/ f\v\^dv = 2N-e^^^ [ [ ffM'dv^dv, 

Ut JRN 4 jRNjf^N 

where 

f 1 - (z/-a) , , 
e^v = / 0{u, a) da = const. 

JsN-i 2 

Notice that, unlike the no-diffusion case, the kinetic energy is not necessarily a 
monotone function of time. However, it is not difficult to show using (jl.lOp (see 
Section 121) that the kinetic energy remains bounded for all times, provided the initial 
distribution function has finite energy. 

Finally, equation ()1.7|) allows us to define the concept of solutions of ()0.1|) which 
we use throughout the paper. Namely, we say that a function / is a weak solution of 
dnUl) if for every T > 0, / G ( [0, T] x R^) , g (/, /) G ( [0, T] x M^) and ^ holds 
for every ip G C^([0, oo), C^(M^)) vanishing for t > T. It can be shown in the usual 
way that if a weak solution is sufficiently smooth (say, continuously differentiable 
with respect to time and twice continuously differentiable with respect to velocity) 
and satisfies suitable decay conditions for large \v\, then it also is a classical solution. 

1.4. Entropy identity. Taking in the weak form ()1.6|) ip = log/ we obtain an 
interesting identity for the entropy Jj^jv f^ogfdv. First, we compute 

Q{fJ)^ogfdv=^[ [ [ f f Jog ^-^\u\b{u, a) da dvdv^ 

' ^Jm.Nj^NJsN-1 J J* 

=11 [ [ ffJlog^-^ + l)\u\biu,a)dadvdv, 
+ ^ / / / {f'fl-ff*)\u\b{u,a)dadvdv*. 

'^JrnJ^n Jsn-1 
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The last term vanishes in the elastic case a = 1; however, as we see below, it is 
generally different from zero if a < 1. To compare the integral of /'/^ to that of //* 
we perform the transformation corresponding to the inverse collision, passing from 
the velocities v', v'^ to their predecessors v and f*. Such a transformation is more 
easily expressed in the variables u and n. Passing to these variables, we can write 
the integral of /'/* as follows: 



1.12) d^-^ I / I ^ J'fl\u-n\dndvdv. 



where S^~^ = {n G S^~^ \ u ■ n > 0}. The "inverse collision" transformation 
{v,v^,,n) {v',vl, —n) has the Jacobian determinant equal to a [10 . Therefore, 
using the first of the equations (jl.ll) . the integral p.l2|) is computed as 

(1.13) d^~^ — I I I f f<t\u ■ n\dndv dv^, 

Changing variables in the angular integral from n to a, we rewrite ()1.12j) as 

(1.14) ~^ I I I ff*\u\b{u,a)dadvdv^ = ^ f f ff^\u\dvdv^. 

« JrN JrN Js^-l a J^N J^N 

In view of fll.llj) and ()1.14|1 the entropy equation becomes 
d 



dt 



f log fdv + 4 



dv 



;i.l5) [ [ ffJlog^-^ + l]\u\biu,a)dadvdv, 

'^JrnJrnJsn-1 V //* //* 

+ ^f^-l) / / ff^\u\dvdv^. 



2 Va2 

In these equations as in all the sequel, the symbol V will stand for the gradient 
operator with respect to velocity variables. Here the first term on the right-hand 
side is nonpositive (notice the inequality logx — x + 1 < 0) and similar to the entropy 
dissipation in the elastic case. The last term in p.l5|l is a nonnegative correction 
term that vanishes in the elastic limit a — > 1. 

1.5. Similarity in the equations and normalization of solutions. As a con- 
sequence of fjl.Sp . the total density (mass) and momentum (mean value) of the 
distribution function are equal to those of the initial distribution. We can write this 
as follows: 

f dv = po = const, and / f Vidv = poVoi = constj, i = 1, . . . , N. 

In fact, we can always assume that po = 1, vq = and p = 1 in ()0.1|) . Indeed, if 
f{v,t) is such a solution to (jO.ll) . then, for every po, vq and p, the function 
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where 

^ = Po^^^f^'^^'^^ and r] = p^^^^^^'^, 
is a solution corresponding to the given values of po; "^o and p. 

1.6. Strong form of the collision operator. Using the weak form ()1.6|1 we can 

derive the usual strong form of the collision operator. We notice the obvious splitting 
into the "gain" and the "loss" terms, 

Q{9j) = Q^{gJ)-Q-{gJ)- 

Assuming that / is regular enough, setting ipiy) = 6{v — vq) in the part of ()1.6|) 
corresponding to Q~{g, /), and using ()1.5|) we find 



Q i9,f)= / fg*\u\b{u,a)dadv^ = f{g*\v\). 

To find the explicit form of Q~^{g, f) we invoke the inverse collision transformation, 
tracing the collision history back from the pair f , to their predecessors, which we 
denote by 'v and 'v^,. Setting ip{v) = 6{v — vo) and arguing similarly to the derivation 
of the entropy identity we obtain 

Q^i9,f)=[ I 'f'g*^\u\b{u,(T)dadv^, 

where '/ = f{'v,t), 'g^ = g{'V:^,t), and the pre-colhsional velocities are defined as 
w 'u w 

(1.16) '^^2"^ 2"' '^*^2"~2"' '^^^^^ 'm = (1 - 7)u + 7 |u|(T, 

and 7 = 

2. Basic Apriori Estimates: Energy and Entropy 

In the classical theory of the elastic Boltzmann equation, the energy conservation 
and the entropy decay are the most fundamental facts which provide the base for 
every analysis. In the present setting naturally we do not have energy conservation, 
and the energy inequality (expressing that collisions do not increase the energy) 
would by no means be sufficient to compensate for that. So the key ingredient will 
be to replace it by the more precise energy dissipation estimate, as follows. 

To study solutions of (jO.lj) and ()().2j) we assume for simplicity that they satisfy 
the normalization conditions of unit mass and zero average; however the estimates 
we derive below will be by no means restricted to such solutions. We use the energy 
equation (jl.iup and apply Jensen's inequality for the last term to get 

3 



f^\u\ dv^ > V— / f{t,v)vdv 
and therefore, 

ff^\ufdv^dv> / f\v\^dv. 
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We then get (in the time- dependent case) the differential inequahty 
(2.1) ^[ f\v\^dv + kif f\vfdv<Ki, 

where Ki = 2N and ki = eN^—r-- Further, by Jensen's inequahty, 



f\vfdv >( [ f\v\^dvy^\ 

V TON / 



and we obtain 

d r „, . , . / /■ „, ,2^ \3/2 

\v\ dv 



dt 



[ f\v\^dv<Ki-ki( [ f 

JRJV V Jrn 

Thus, if f\v\'^dv > {Ki/kiY^^, then ^ J^^^ /|f prff < and so, 
sup / / |f p c/t> < max I / fo\v\^ dv, [Ki/kiY^^\ . 

t>0 JrN l J^N J 

In the steady case the derivative term drops in (j2.ip . and we obtain 

Let us introduce the following weighted spaces: 
(2.2) L^(M^) = {/ I f{v)' e Li(M^)}, 

where > and {v) = (l + |fp)^/^. We then define the norms in as Jj^jv \f\{v)''dv, 
which for / nonnegative coincide with the moments J^f^ f (v)'' dv. The above ar- 
gument implies apriori estimates for the steady solutions in L3(M^), and for the 
time-dependent ones in L°°([0, oo), L^(M^)) and Lj„^([0, oo), L^(M^)). We empha- 
size that the bounds depend on a and deteriorate in the elastic limit a —>■ 1. In fact, 
these bounds for a < 1 make a most striking contrast with the classical Boltzmann 
equation for elastic particles. 

Next, using the entropy equation ()1.15|1 we show that the entropy is bounded 
uniformly in time, for initial data with finite mass, kinetic energy and entropy. To 
obtain this, we first estimate the second term in ()1.15jl using the Sobolev embedding 
inequality: assuming for simplicity here that > 3, we have 

ivv7rrf^>cii/iL.*, 



where p* = N/{N — 2). Further, we have the inequality 

(2.3) / flog fdv< Cell fill,*, 

jRN 

for all £ > 0. Indeed, obviously, for every 6 > 0, 

(2.4) [ f log fdv<Cs[ f^'dv. 



Further, by Holder's inequality, for 5 < p* 
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where v = , , ^,1 , . Therefore, 

(p*-l)(l+<5) ' 

which together with ()2.4|) imphes ()2.3|) . Now, coming back to estimating the terms 
in the entropy equation we get 

(2.5) ^/ f\ogfdv + cJ[ f\ogfdvy^'<C\\f\\l.. 

The estabhshed bound in L°°([0, oo), L2(M^)) imphes that for initial data with finite 
mass and energy, the right-hand side of ()2.5j) is bounded by a constant, and we obtain 
by Gronwall's lemma, 

sup / f log f dv <C{J^^ fo\ogfodv,\\fo\\Li) . 

Integrating ()1.15|1 in time, we also get v7 ^ L'^il'^^T], H^(R'^)), for every T > 
0, which implies in particular, / G L^*([0,T] x M^), where the constants in the 
estimates depend on the initial mass, energy and entropy of the solutions. For the 
steady solutions we obtain a particularly simple estimate ||Vv7||l2 < C||/||/^i. As 
the reader will easily check, our assumption that > 3 is just for convenience, and 
can easily be circumvented in dimension 2 by the Moser-Trudinger inequality, or just 
the local control of all norms of / by ||Vv7||l2, together with a moment-based 
localization argument. 

3. Moment inequalities 

We further look for apriori estimates of the solutions in the spaces Ll ()2.2j) with 
k > 2. Such estimates will play a very important role in our study of regularity, 
which we perform in Section |3 The key technique for obtaining the necessary 
estimates is the so-called Povzner inequalities jSH [13 1121 IHH HI IHOj which we here 
extend to the inelastic case. 

3.1. The Povzner-type inequalities. We take V^(a;), x > to be a convex non- 
decreasing function and look for estimates of the expressions 

(3.1) q mv, v., a) = ^{\vf) + ^(|<p) - ^(1^;^) - ^{\v^\') 
and 

(3.2) q[4;]{v,v,)= [ (^(\vf)+^(\v'f)-i;{\v\')-ij{\v,\'))b{u,a)da, 

which appear in the weak form of the collision operator p.6|) . 
Our aim is to treat the cases of 



(3.3) 



ipi^x) = x^, and iIj{x) = (1 + xY — 1, p > 1, 
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and also truncated versions of such functions which will be required in the rigorous 
analysis of moments in Section El Thus, we will require functions ip to satisfy the 
following list of conditions: 

(3.4) -^{x) > 0, X > 0; V^(0) = 0; 

(3.5) ip{x) is convex, C^([0, cxd)), ip'^x) is locally bounded; 

(3.6) ip'iax) < r]i{a) -^'{x), x > 0, a > 1; 

(3.7) ilj"{ax) < r/2(a) ^p"{x), x>0 a > 1, 

where rii{a) and 772(0) are functions of a only, bounded on every finite interval of 
a > 0. The above conditions are easily verified for the functions ()3.3p . 

We will further establish the following elementary lemma. 

Lemma 3.1. Assume that ip{x) satisfies (j3.4j) - (j3.7j) . Then 

(3.8) ipix + y) - tpix) - ij{y) < A{xilj'{y) + yij'{x) ) 
and 

(3.9) ipi^ + y) ~ i^ix) — ipiy) > bxyip"{x + y). 
where A = 771 (2) and b = (2772(2))"^ 

Proof. To establish the first of the bounds assume that x > y. Then, since ip{y) > 0, 
ipix + y) — '>p{x) — 'tpiy) < ip{x + y) — ip{x) 

ry 

ip'ix + t)dt < / r]i{2)ilj'{x)dt = Ay^'{x), 



By symmetry we have 

i^{x + y) -^p{x) -i^iy) < Axip'iy), 

when X < y. This proves the required inequality for all x and y. To prove the 
second of the bounds in the lemma, we can write, using (|3.7p and the normalization 
^(0) = 0, 

ry ry rx 

ip{x + y)-ij{x)-'ilj{y)= / {ij'{x + t) -ij\t))dt = / / ij"{t + T)drdt 



^0 



y rx 



>{V2{2))''^"{x + y) / / X{t+r>i.+y)/2}dTdt={27]2{2)y'xynx + y). 

Jo Jo 

This completes the proof. □ 

In the sequel, we shall use some relations involving post-collisional velocities v' and 
v^. It becomes more convenient to parametrize them in the center of mass-relative 
velocity variables. We therefore set 

/QinN / w + AImIu; , , w — X\u\uj 
[o.iO) V = and = , 
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where w; = u + f*,M = f — f*, and u is a parameter vector on the sphere S'^~^ (see 
Figure 1). We have 

where P = and u = u/\u\, and therefore, 

(3.11) A = A(cosx) = (1 COSX+ V(l -/3)2(cos2x- 1) 

where x is the angle between u and u. Notice that 

< a < A(cosx) < 1, 

for all X- With this parametrization we have 

|wp + A^lwp + 2A|m| cos/x 



4 



(3.12) , ,„ „, , , 

I /|2 |w| + A — 2A|'u| |w| cos/i 
\vj = ^ , 

where is the angle between the vectors w = v + v^, and u. 

Lemma 3.2. Assume that the function ip satisfies ()3.4j) - ()3.7j) . Then we have 

q[ijj] = -n [ijj] 

where 

p[iP]<A{\v\^^'i\v4') + \v,\'^'{\v\^)) 

and 

n [^] > k(A, fi) {\v\^ + \v,\^)^ \v\^ + \v,\^). 
Here A is the constant in estimate (j3.8j) . 

«^(-^>Ai) = ^-^^ ('72(A"^))"^ sinV, 
and b is the constant in estimate (13. 9j) . 

Proof. We start by setting 
and 

The estimate for p [ip] follows easily by ()3.8|) . It remains to verify the lower bound 
for n [ip] . For this we use (j3.9|) , noticing that ijj is monotone and that 1 p + | f p > 
l"^'!^ + l"^*!^- We then obtain: 

nm>H\vf + \v'f)-ij{\vf)-iji\v'f) 
>b\vf\v'fr{\vf + \vT) 



bav',v:)i\vf + \vXfr{\vf + \v'X) 



where 
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Further, using ()3.12p . we get 

1 / 4A ^ w 



4 

Finally, noticing that 



4A2 




2 




2 




■u 


2 + 




2^ 



we obtain 



I /|2 , I /|2 A2|m|2+|w;|2 ,2l^P + kP x2 n |2 , 1 |2n 



nm>^ X' (r/2(A-'))-^ sinV + \v^\y f'(\v\^ + \v^\^). 



This completes the proof of the lemma. □ 

Lemma 13.21 gives us the basic formulation of the Povzner inequality for the con- 
sidered class of test functions In the example ipi^x) = we have p [ip] ~ 
C(|t;p|f*pP"^-|- pP"^) and n[ip] ~ c(|wpP+ outside the set where k(A, /i) 

is small (which amounts to a small set of angles). This implies that that the nonpos- 
itive term —n [ip] is dominating, at least when |t>| >> |t>*| or |t>*| >> |t>|, which are 
the most important regions of integration from the point of view of calculation of 
moments (cf. also fT^ I31j). We can further simplify the inequalities and get rid of 
the dependence on the angular variables, by integration with respect to cr G 5*^^^. 
We then obtain the following lemma. 

Lemma 3.3. Assume that the function if) satisfies ()3.4|) - fl3.7|) . Then 

q[ij] < -A;(|^;|2 + |^;,|2)2^"(|^;|2 + |^^|2)+ ^(|^|2^'(|^^|2) + |^^|2^'(|^|2))_ 

where the constant A is as in Lemma 1,9. M and k > is a constant that depends on 
the function ip but not on a. 

Proof. For the proof we notice that A(cosx) is pointwise decreasing as a \ and 
so, 

A(cosx) > cosx, for cosx > 0, 

for all a > 0. We then denote cos 6* = {u ■ a), bo{cos9) = b{u, a), and estimate the 
integral 

(3.13) / k{X, fi)bQ{cos9) da > / k(A, /i) 6o(cos 6*) (icr, 

J gN-i J{cosx>£0: sin/^>ei, 1— cos6>£2} 

setting ^1 and 62 small enough. The integrand on the right-hand side of ()3.13|) is 
bounded below by a constant, and so is the area of the domain of integration. (The 
verification of the last statement for the condition sin /i > ei is somewhat tedious 
and is achieved by changing the variables of integration from uj to a: we omit the 
technical details). We therefore find that the integral ()3.13|) is bounded below by a 
constant k > 0, independent on a. The rest of the claim is easy to verify. □ 
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Finally, we present estimates for the integral expression ()3.2|) multiplied by the 
relative speed, in the cases when ip{x) is given by one of the functions ()3.3|) . 

Lemma 3.4. Take p > 1 and iplx) = x^. Then 

\u\q[iP]{v,v,) < -kp{\v\^P+' + \v4^''+') + Ap{\v\\v,\''^ + \v\^''\v4). 
Also, take iIj{x) = (1 + x)^ — 1, then 

\u\qmv,v,) < -k,{{v)'P+' + {v,)'^+') + A,{{v){v,fP + {v)'p{v,)). 
Here the constants hp and Ap are independent on the restitution coefficient a. 

Proof. We use Lemma 13.31 and the inequalities 

\\v\ — \v^\\ < \u\ = \v — < \v\ + 
Then in the case iplx) = x^ the bounds have the form 

-p{p - l)kp \u\ + Iv^iy + pAp \u\ + \v\''''"^\v,f) • 

The terms appearing with the negative sign are estimated using the inequality 

\u\ i\v\^ + \v,\^Y > ^i\v\^^'^' + Iv]^"^') - ^Mv*]^" + ■ 

For the remaining terms we have 

\u\ {\v\^\v,\^P-^ + < Cpi\v\\v,\^P + Ivl^^M , 

which completes the proof of the first part of the lemma. The case ijj{x) = — 1 

can be treated by arguing along the same lines, by using the inequalities 

{\v\' + \v.\'f>^{l + \v\' + \v,\y-l 

and > (1 + |f - 1. □ 

3.2. Estimates for higher-order moments. The Povzner-type inequalities of 
Lemma EiH allow us to study the topics of propagation and appearance of moments. 
We find that results known for the classical Boltzmann equation with "hard-forces" 
interactions ^3 ^] transfer to present case. We introduce the notation 



Ut)= / fivydv, 

and denote by Yg the corresponding steady moment. 

Lemma 3.5. Let f be a sufficiently regular and rapidly decaying solution of (lU.ip . 
Then, the following differential inequality holds: 

(3.14) + 2A;,r,+i < Ks{Y, + ^.2) 

dt 

where Kg and kg are positive constants. Further, 

snpYsit) < Y; = max{y;(0), (iT./fc,)'}, 

i>0 
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and for every r > 

(3.15) fYs+i{t)dt < ^-^ + ^/V ;. 

Jo l^s 

Finally, for the steady equation (j(J.2j) we obtain the apriori estimate 

Proof of Lemma \3.f\ Using the weak form of equation ()n.l|) with ip{v) = (v)^ we find 

(3.16) 4/ fi^ydv- [ Af{vrdv= [ Q{fJ){vydv. 



dt 

Estimating the moments of the colhsion integral according to Lemma f3. 41 we get 



(3.17) 



Q(/, /) {vY dv < -2ksYs+i + 2AsYiYs. 
The moments of the Laplacian term are computed as follows: 

fA{vYdv= [ f{{s{s-2) + sN){vY-^-s{s-2){vy-'')dv 

= {s{s - 2) + sN) - s{s - 2) F,_4. 

Combining ()3.16|) and ()3.17|) and neglecting the non-positive Ys_4 term, we obtain 
inequality (pTTijl with Ks = max{2As, s{s - 2) + sN}. 

To obtain a uniform bound for Ys{t), we use Jensen's inequality to write 

Ys+i > (n)(^+^)/^ 

Then we find, estimating the right-hand side of ()3.14j) by 2KsYs, 

^n<-2A;,(n)(^+i)/^ + 2ir,F,. 
dt 

Thus, Xj'(t) < if > {Ks/ksY, and so, the upper bound for sup^>QFs(i) niust 
hold. 

Further, integrating in time we obtain 

2k, I < 2k,ty: - Y{s) + r(o) < {2K,T + i)y;. 



which proves ()3.15|1 . 

Finally, the last inequality is obtained by the same arguments as ()3.14|1 applied to 
the steady equation. □ 

Based on the Lemma just proven we can make the following conclusions about 
the behavior of the moments of the solutions. First, if a moment Y^ is finite initially, 
it propagates, that is, it remains bounded for the whole time-evolution. Further, the 
integral condition on Kj+i implies the appearance of moments of order s + 1: these 
moments become finite after arbitrarily short time, even if they are initially infinite 
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(cf. [12] )• Indeed, suppose that Fs+i(0) = +00, then for every r > there is a 
to < T such that K;+i(to) < +00. Then, applying the Lemma to Kj+i, starting with 
t = to, we obtain that for every to > 0, 

supK,+i(t) < Cto,s, 

t>to 

which imphes the above statement. The last part of the Lemma implies an important 
statement concerning the moments of the steady solution: on the formal level, every 
solution that has a finite moment of order s > 2 has finite moments of all positive 
orders. In fact, in view of the Ll{R^) estimate of the previous section, this implies 
that every solution with finite mass should have this property. 

4. LP BOUNDS AND APRIORI REGULARITY 

In this section we study the apriori regularity of solutions to (jO.lll and ()().2|1 . The 
presence of the diffusion term in the equation makes it plausible that solutions to 
the steady equation should be smooth, and those for the time-dependent equation 
should gain smoothness after arbitrarily short time. However, to realize this idea 
we need to make use of the particular structure of the collision term. As we will see 
below, the moment bounds of the previous section will also be of crucial importance. 
We start by establishing the bounds for the collision operator in the spaces with 
a polynomial weight, extending the results well-known in the case of the classical 
Boltzmann equation, and first derived by Gustafsson [21]. Below, we shall establish 
these bounds by adapting the simple strategy that was suggested in jHHl Chapter 2, 
Section 3.3] and later developed in to establish improved bounds in the elastic 
case. 

4.1. LP bounds for the collision operator. We will use the following weighted 
LP spaces: 

L^(M^) = {/ I f{v)' e L^(M^)}, 

where (v) = (1 + |f p)^/^. The necessity to introduce a weight comes from the 
presence of the factor |m| in the hard sphere collision term p.6p . The collision 
operator is generally unbounded on LP; in order to control its norm we will invoke 
the L^ norms with higher powers of (v). The precise formulation of this statement 
is given in next lemma. 

Lemma 4.1. For every 1 < p < 00 and every k > 0, 
where C is a constant depending on p, k and N only. 
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Proof. We fix an exponent 1 < p < oo. It is easy to estimate the "loss" part 
Q'i9,f) = i9*\v\)f, using the inequality 

\9*\v\\< II^IIli(^), 

from which it follows 

(4.1) ||Q-(^7,/)IIl^<II^7||liII/1Il;:,,- 

We now turn to estimate the term: starting from the weak form we find 

\\Q^{9j){v)'\\Lr'= sup / Q+{gJ)^{v)'dv. 

(4.2) 

= sup / / f 9*\u\ ip' {v')'' b{u, a) da dv dV:f,. 
\m^p,=iJm Jm is'^-i 

By using the inequalities \u\ < (f ) + (f*) and {v')^ < {{v) + {v^))'^ the integral ()4.2|) 
is bounded as 

(4.3) / / f9*i{v) + {v*)f^^ [ ^'h{u,a)dadvdv,. 

We now see that the problem comes down to estimating the integral 

S[il!]{v,v^) = / ip'h{u,a)da 



in either L~(R^, L?''(M^)) or L'^iM.^^^, L^' {M.^)). In fact, we split S[^] into two 
parts 5"+ [^] and S- [^] and prove the bounds for each of the parts in the respective 
spaces. We set 



S±[i'\{v,v^) = I ilj'b{u,a)da 

J{±ua>0} 

and establish the bounds for 5*+ and S- in the following proposition. 

Proposition 4.2. The operators 

S+ : L'?(M^) ^ L°°(Mf ,L'?(M^)), 

are bounded for every 1 < q < oo. 

Proof. We prove the L'^ bounds by interpolation between L°° and L^. The L°° 
estimates are clear due to the boundedness of the domain of integration. To check 
the bounds we assume without loss of generality that ^Z' > and calculate the 
norms as follows: 

\\S-[iP]{v,v^)\\l^kn-)= [ [ 'ip(v + ^{-u+\u\a))b{u,a)dadu 

+ z) [ H^(^^ ^U)X{(...(..))<o> 
JsN-i \J^{u{z,a),a)\ 
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Here, z = v' — v = | {—u + \u\a) , and J_(m, a) is the Jacobian of the transformation 
u ^ z (for fixed a): 



2) V \u\ 

The condition -u ■ cr < ensures that | J_| is bounded below by (f and then, 



\\S^[^Ij]{v,v^)\\l^^nj < (^2 7 j^^_h{u,a)da\\i)\\Li ^ Vi) 

for every f G M^. 

Similarly, for the 5*+ term we have 

\\S+[i)]{v,v^)\\L^(^M.N)= I I ^Ij (v^ + ]-{{2- P)u + p\u\a))h{u,a)dadu 

Riv y^iv-i |J+(u(2;,o-),cr)| 

where now z = f ' — f.,, = ^((2 — I3))u + /3 |M|cr), and 

jAu,a) = {—^] 1+ ^ 



2 / V 2-/? |m| 
Then, since {u ■ a) > 0, we can argue similarly to the previous case to obtain 



\S+[i)]{v,V^)\\L^^N) < 



2-/3 



-N 



uniformly in G M^. The statement of the proposition now follows by the 
Marcinkiewicz interpolation theorem. □ 

End of proof of Lemma \41\ Combining the bound (j4.3|) with the ones proven in 
Proposition 14.21 we find 

< I I f9.{{v) + {v:)f^^[S+[ij]{v,v,) + S_[ij]{v,v,))dvdv, 



<Ck g. f{{v)''+' + {v,)''^')S4^]{v,v,)dvdv, 
+ Ck[ f[ g.{{v)'''-' + {v,)'+')S.mv,v,)dv,dv 

< c ( WghAlfhi^, + MliJIMl. + llfhAlghi^, + ll/IL^JklU.). 

since HV^H^.p' = 1- From this the conclusion of the lemma follows easily. □ 



ON THE BOLTZMANN EQUATION FOR DIFFUSIVELY EXCITED GRANULAR MEDIA 19 

4.2. regularity: Steady state equation. We start by establishing apriori 
estimates for solutions to the steady equation ()0.2|) . for which the analysis is per- 
formed in a rather more direct way than for the time-dependent problem. We first 
show the bounds in the Sobolev spaces with the weight (v)'' = (1 + |f p)'^/^: 

The main tools are the coercivity of the diffusion part, the estimates of the collision 
operator in L^, and the interpolation inequalities for spaces. The constants in 
the estimates are expressed in terms of the moments. In all this section, we 
shall assume for simplicity that > 3, but there is no difficulty to adapt the proofs 
to cover the case = 2 as well. We begin with an estimate for the gradient in L^. 

Lemma 4.3. Assume that the function f G H^iM^) fl Lj(]R^), where r = is 
a solution of ()().2|1 . Then 

Ilv/|U2<CA5^ 

where 



A=\\f\\Li, B=iijui, 
and C is a constant depending on the dimension. 

Proof. Multiplying equation ()0.2j) by /, integrating and applying Holder's inequality 
yields 

(4.4) [ \Vf\'dv= [ Q{fJ)fdv<C\\fh4Q{fJ)\\ 



LP 5 

for all 1 < p < oo. We choose p = 2* = 2N / {N — 2), where 2* is the critical Sobolev 
exponent, and apply the Sobolev's embedding inequality 

(4.5) ll/IL- <C||V/|U^ 

(Note: for = 3, 2* = 6 and (2*)' = 6/5.) Then, by Lemma EU 

(4.6) ||Q(/,/)IL,.*)' <C||/||^(..y||/|U;. 

To estimate 11/ 11 (2*)' we use the following interpolation inequality for weighted 

^1 

norms {(p is any weight function), which can be easily verified using Holder's 
inequality: 

(4.7) \\f^%.<\\f^'TL^A\fv''\\l-..', 

where 

1 = - and kiu + A;2(l — z/) = k. 

Qi 92 q 

Now, interpolating the norm in for q = (2*)' between qi = 2* and q2 = 1, we get 

(4.8) ii/ii,(.v< 11/11^.* ii/iiir 
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where v and r are determined from the following equations: 

V \ — V 1 1/1 \ 1 

1 = - — ^ and r (1 — z/j = 1, 

2* 1 (2*)' ^ ' 

so that 

, , , 1 A^ + 2 

(4.9) ^=ivT^ ^ = T3^ = ^- 

Combining estimates ()4.4|) - ()4.8p we obtain the inequality 

(4.10) ||V/||i. < C7||/|U.||/||i7l|/||itr < C5Ai-1|V/||ir, 

from which the conclusion of the lemma follows. □ 



The result of the Lemma implies a bound for the solutions in the space H 
Indeed, by the Sobolev embedding, 



1^2. < CAW. 

Interpolating between V- and L^* using inequahty ()4.7|) we get a bound for the 
I? norm, which then implies a bound in . Since the constants in the estimates 
depend on the L\ norms only, and the latter are controlled by the moments bounds, 
we gain an apriori control of the norm by means of the mass and the energy 
only. We next see that the derivatives of the solutions have an appropriate decay, 
so even L\ norms for all A; > are bounded. 



Lemma 4.4. Let j he a solution of equation ()0.2|) and assume that f G Hl(R^) fl 
L[^^y^,(^^); ^here k>0 andr = ^. Then 

\\V{f{v)')h2<C{A^A; + e^A,), 

where 



^1 = 11/11 W ^2 = II/II^U.' ^3=||J||Li_,^^. 

and C is a constant depending on the dimension N. 
Proof. Integrating equation ()().2j) against f {v)'^^ we obtain 
(4.11) / Vf-VU{vY')dv= I Q{f,f)f{v)''dv. 



Using estimates from the previous lemma, the right-hand side can be bounded above 
as follows: 

IIQ(/,/)(^)1L(^*,'II/(^)1l- 

(4.12) 

^ ^11 J" /..\ fc+i II II J- /..\ fc+i II ll_f/,.\fcl, 



<C\\f{v)''^Xi-^.y\\f{v}'^^%4f{vy 



Interpolating as in ()4.7p we find 



(4.13) ll/(^)'+^llL(-)'<ll/III-ll/lliT'^ , 
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where u and r are as defined in ()4.9|) . Therefore, combining ()4.12|) with ()4.13|) we 
bound the right hand side of ()4.11|) by 

cdl/lliT'^ ll/llL^,JII/(^)1litr 

<CA\-^A,\\WU{v)'Wj^.r 
The integral on the left-hand side of ()4.1H) is estimated as follows: 

/ Vf-V{f{vf'')dv= [ \V{f{v)'')\^dv- [ f\V{v)''\^dv 



(4.15) 



>l|V(/(t')'^)||i.-A;^||/||i. ^ 

fc — 1 

Further, interpolating the Ll_-^^ norm between L^* and we get 

\\f{v)'''%. < c\\f{vm.4f{v)'%-.' < c\\v{f{vmu\f\\i.\ 

'=2 

where 

A 1 - A 1 , , N l + v 

(4.16) \ = -, so that A = — = , 

^ ' 2* 1 2' A^ + 2 2 ' 

and 

A; - 1 = Afc + (1 - A)A;2, so that k2 = k-^ — - = k-2/r. 
Gathering the above inequalities and noticing that 2A = 1 + z/ we obtain: 



\Vif{v)')\\l<C{Al-'^A, + eAl~'^)\\V{f{v)'''^^'' 



Dividing by the norm of the gradient to the power 1 + u we get 

\\V{f{v)')\\L2 < {CAl~''A, + k'Al-^)^. 

Noticing that = r and using the inequality (x + yY < Cr{x^ + y^) we arrive at 
the conclusion of the lemma. □ 

Using the Lemma just proven we find bounds for solutions / in for every 
k > 0. Indeed, using the inequality 

K'^niv)']' <C{\V{f{v)')\' + \fV{v)'^\') 

and interpolating in the second term between L^* and Ll_2/r we get 

iiv/iii^ < c{\\v{f{vf)\\i. + iMfivmmni-,'^ ), 

^ k — 2/r 

from which an estimate in terms of the moments follows. Further, by interpola- 
tion inequality ()4.7p . 

A II r||l-A 



-^fe/(l-A) 

and so, in view of our earlier remarks, the norm in is also estimated in terms of 
moments only. Summarizing the results obtained so far, the solutions are controlled 
apriori in if^(R^) for any > in terms of mass and kinetic energy only. 
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4.3. Schwartz class regularity: Steady problem. Our next aim is now to es- 
tablish a priori bounds for solutions to ()0.2j) in the spaces 

H-iR'') = {fe I V™/ G L^(M^), 1 < m < n}, 

for all 1 < < CX3 and all < A; < oo. We use induction on n, differentiating the 
equation in v in each step. The base of the induction is given by Lemma 14.41 We 
recall the following rule for differentiating the collision integral. 

Proposition 4.5. Let f and g be smooth, rapidly decaying functions of v. Then 

'^Qi9J) = Qi^9J) + Q{9,'^f)- 

Proof. We use the splitting into the "gain" and "loss" terms, Q{g, f) = Q~^ig, f) — 
Q~{9, /). Since Q~{g, f) = f{g * \v\), the differentiation rule for the "loss" term is 
obvious. To prove the proposition for the "gain" term Q~^{g,f) we represent it as 
follows, using 

Q^i9J)=f [ f(v+ ^ ^ ^ )g(v- ^^^ ) -^\u\b{u,a)dadu. 

Since 'u is a function of u and a only, the statement follows by differentiation under 
the integral sign. □ 

Remark. The above statement is in fact a corollary of the following abstract state- 
ment which can be proven very easily: Let Q be a bilinear operator commuting with 
translations, continuously differentiable; then VQ{g,f) = Q(Vf,g) + Q{f,^g)- 
Thus, the differentiation formula of Proposition 14.51 can be seen as a consequence of 
the translation invariance of Q. 

As a direct corollary of Proposition 14.51 higher-order derivatives of Q can be 
calculated using the following Leibniz formula: 



d'Q{gJ)= (])Qi9'-'9,d'f), 



0<l<j 

where j and / are multi-indices j = (ji . . . jjy), and / = {li . . .Ijy); 

u u^_^ . . . u^^ , 

and (-J) are the multinomial coefficients. Thus, for every multi- index j, by formal 
differentiation of ()0.2j) we obtain the following equations for higher-order derivatives: 

(4.17) -A9V= {^)jQid'~'9,d'f)- 

o<i<j ^ ^ 

By applying the methods developed in Lemmas 14.31 and 14.41 to equation ()4.17p we 
arrive at the following result. 
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Lemma 4.6. Let f he a solution to (jHS)), such that f G HJ!:^^(R^), with n > 0, 
k > and > 1 + y • Then 

||v"+VIIli <c{i + k + i) (1 + Wv-fhi J, 

where C is a constant depending on n and N only. 

Proof. Taking a multi- index j with |j| = n, multiplying equation ()4.17|) by f {v)'^^ 
and integrating by parts we obtain: 



(4.18) 



Vd^f ■ V(9V {vf^) dv 
til-. V/ Jm 



0<l<j 

Similarly to ()4.15|1 . the left-hand side can be written as 

\Md^f{vn\h-\\d^fV{v)''\\h. 

Each integral on the right-hand side of ()4.18|) can be bounded above by using 
Cauchy-Schwartz's inequality and Lemma f4. II as follows: 



g(9^-7, d^f)&f{vf^ dv 

< md^^'f.d'f)\\Ll\\d^f\\Ll < c||9'/|Li^j|a'-vilL|,J|vvili.^J|5VILi. 

Now, the norms can be estimated as follows: 

Pf\\LU,<m'-'\\L^ \\d'f\\LU^<C\Wf\\,.^^ 

as soon as /i > 1 + y. Gathering the above estimates we obtain: 



2 

L2 



< 



wfy{vni.+c\\cPf\\,. y: (])\\d'f\\Lijd'-'fu,, 



0<l<j 



<e\\d^f\\l, +c{i + \\f\\i„..wf\\l. +c{i + \\d^f\\i. 

Since 
we obtain 

II v^viii. <c{i + e + (1 + ii^viii. ). 

Taking the sum over all j with \j\ = n implies the estimate of the lemma. □ 

Lemma 14.61 gives us a way to estimate higher-order derivatives of solutions in 
terms of lower-order ones. Thus, provided we have a solution to ()0.2|) that has all 
Hi norms bounded in terms of mass and energy (as we assumed in the previous 
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section), we can derive bounds in if| for every k, and then proceed by induction, 
obtaining bounds in if^(M^), for all n and all A; > 0. We then obtain 

n>l,fc>0 

where S is the Schwartz class of rapidly decaying smooth functions. Notice that the 
bounds in each of the spaces if^(]R^) can be expressed in terms of mass and energy 
of the solutions. 



4.4. Regularity for the time-dependent problem. An analysis of the regularity 
of the time-dependent solutions can be performed in the same vein as for the steady 
problem. Using the estimates obtained in the previous section in combination with 
Gronwall lemma will give us results for the time-dependent equation (lO.lj) . Our first 
lemma is an analog of Lemma 14. HI 

Lemma 4.7. Let f be a sufficiently regular solution to (lO.lj) with the initial condition 
/(■,0) = /o G L^(R^), such that f has a moment of order r = bounded 
uniformly in time. Then 

\\f{;t)\\L2<C, 0<t<cx), 

and 

for every < T < oo, where the constants C and Ct depend only on N , ||/o||l2 and 
sup||/(-,t)lUi.. 

Proof. Integrating equation (j(J.l|) against / we get, arguing similarly to the case of 
the steady problem: 

(4.19) ^|l|/|li. + l|V/||i.<ir(t)||V/||ir, 

where K{t) = C \\f\\i^i \\f\\]^'' and u = as in ()4.9|1 . By interpolation and Sobolev 
embedding, 

ll/llL2<ll/lli-ll/lliT'<C^I|V/||i.||/||i7^ 
where A = as given by ()4.16|) . Therefore, 

(4.20) ||V/||i. > k\\f\\%\ 

where k = is a constant. Distributing the term ||V/||l2 in (|4.19p 

equally between the left and the right-hand sides and using inequality ()4.20p we 
obtain 



(4.21) + k\\f\\%' < -||V/||i. + 2ir(t)||V/||ir 



ON THE BOLTZMANN EQUATION FOR DIFFUSIVELY EXCITED GRANULAR MEDIA 25 

The function X i— — + 2K{t)X^^'^ , appearing on the right-hand side of ()4.21|) 
has a global maximum (1 + z/)^''^^(l — v)K{t)'^'^ = CK{ty^ , so we obtain 

(4.22) + k\\f\\%' < CK{tr < CK'\ 

where K = sup-ft'(t) < sup Applying a Gronwall's lemma argument to (j4.22|) 

t>o t>o 

we then obtain a bound of the norm of / in terms of ||/o||l2 and sup||/||2,i. 

t>o 

Further, integrating ()4.19|) over time, we get 

which proves the second claim of the lemma. □ 

Similar results can be established about the time-dependence of the norms of 
the solutions. 

Lemma 4.8. Let f be a sufficiently regular solution to with initial data fo € 

Ll{M.^), where k > 0, and such that f has a moment of order r{k + 1), where 
r = ^^Y^, bounded uniformly in time. Then 



||/(-,t)|U2 <C, 0<t<oo, 

and 

liv/ll 

for every < T < oo, where the constants C and Ct depend on N, ||/o||l2 and 
sup ||/(-,t)||Li only. 

Proof. Multiplying the equation by f{v)'^^ and integrating we obtain: 

^^11^11'?+/ ^f-'^if{vr)dv= [ Q{fj)f{vrdv 

Following the steps of the proof of Lemma and distributing the term || V(/(f )'^) |p 
evenly between the left and right-hand sides we obtain the following differential 
inequality: 

(4^23) |ii^<^^'ii'^ + w'^ifivnih < -\\^if{vn\h 

+ (cA^ity-'^A^it) + eA.ity-'') \\vif{v)'^)\\ir, 

where Ai{t), A2{t), and A^lt) are the moments defined in Lemma f4. 41 The uniform 
bounds of the moments imply that the right-hand side of (j4.23|) is bounded above 
by a constant. The left-hand side is estimated below as 

d 

dt' 



■.\\f{vnh+c\\f{vnf 
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analogously to ()4.19|) . Thus, by a Gronwall-type argument we obtain that the L^- 
norm of / is bounded uniformly in time. Integrating ()4.23p over time we also get 
the second claim of the lemma. □ 

Finally, we establish the following analog of Lemma 14. (jl which will allow us to 
study the regularity of higher-order derivatives. 

Lemma 4.9. Let f be a solution to (jO.lll with initial data fo G HJ!(R^) where k > 
and n > 0, such that f has a moment of order r* = r{2"'{k + yu) — 2yU + 1), where 
r = and /i > bounded uniformly in time. Then 

and 

II/IIl2([0,T],H^+1(R'V)) < Ct, 

for every < T < oo, where the constants C and Ct depend on N , \\fo\\Hi o-nd 

sup only. 

t>o 

Proof. We will use induction on n. The case n = is already proven in Lemma 
14.81 Assuming that the statement of the lemma holds for n — 1, we differentiate 
the equation in v and argue as in the proof of Lemma 14. 6[ obtaining the following 
inequality: 

^;^liv"/lli. + Iiv"+Vlli. < c{i + e + \\f\\l.-.){l + ||V"/llk ) 

We estimate || V"/||^2 integrating by parts and using Young's inequality (cf. [T^): 

liv^/lli. < 5||V"+Vlli. + Cs\\V-'fhi . 

Then, since we assumed / to be bounded in H2(k+iJ.) 

V"/lli| < -II V"+Vlli. + C,{n, fc)5|| V"+Vlli. + C,{n, k, 6). 
Choosing S suitably small, we obtain the conclusion by Gronwall's lemma. □ 

Lemmas 14 . 71 - allow us to make the following conclusions about the regularity 
of solutions to (jO.ll) . Provided a sufficient number of moments is initially available, 
the regularity of the initial data is preserved with time. Moreover, the established 
bounds for the derivatives in L^([0,T] x M^) imply that after arbitrarily short time 
the derivatives d^f{-,t) of any order are in L^(]R^), and then they propagate in 
time. Thus, on the level of apriori estimates we find that the solutions become 
immediately infinitely smooth in v and decay faster than any negative power for |f | 
large. 

We can also see that the solutions are infinitely differentiable in t. Indeed, in view 
of the established HJ^ regularity we have /(■, t) G 5(]R^) for t > 0, and then equation 
()0.1|] implies dtf{-,t) G 5(]R^), for every t > 0. Differentiating the equation in 
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time and proceeding by induction we find also that (9™/(-,t) G 5(]R^), for every 
m = 1,2,..., and for every t > 0. The time derivatives also remain bounded 
uniformly in time. 

5. Existence 

We next proceed with a rigorous proof of existence that will also justify the formal 
manipulations performed in the derivation of apriori inequalities. 

Theorem 5.1. For every /o > 0, /o G fl LlogL(R^) there exists a nonnegative 
weak solution 

f G L-([0,oo), Li(M^)), /log / G L-([0, oo), Li(M^)) 

to equation fl0.1|) . with the initial condition f{-,0) = /q. Furthermore, if in addition 
fo e Lin L2(M^), where r = max{2, then for every to > 0, 

/GCr([to,oo),5(M^)), 

where denotes the class of functions with bounded derivatives of any order, and 
S is the Schwartz class of rapidly decaying smooth functions. In particular, for 
t > 0, f is a classical solution of (jO.lj) . 

Theorem 5.2. For every p > there exists a nonnegative solution f to ()0.2|) . 

/ G 5(]R^), satisfying / f dv = p. 

Jm 

Furthermore, every nonnegative solution in L]. fl L^(M^)), where r = max{2, ^^^} 
is in fact in S(R^). 

Proof of Theorem \5.1\ We assume that the initial datum /q is in C°°(]R^) and has 
compact support (we will remove this assumption in the end of the proof). We also 
introduce a truncation in the collision term by replacing the factor |m| in by 

(5.1) lulm.A^ = + M} 

and m > 0, M > are truncation parameters. We then denote by Qm,ivi{,f, f) the 
corresponding collision operator. 

The first step of the proof will be to find approximating solutions which we define 
using the following truncated problem 

9J-A,/ = Q^,m(/,/), ^;GM^ tG[0,T], 
f{Q,v) = fo{v), 

where m, M and T are fixed positive parameters. We will denote by / solutions to 
()5.2|) . keeping in mind that they generally depend on m and M. 
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The solutions will be constructed by applying a fixed point argument to the fol- 
lowing approximation scheme: 

dtf-AJ + Mf = Q,n,Mi9,9) + Mg, veR', t e [0,T], 
f{v,0) = fo{v). 

Here g is a nonneg ative function from L°°([0, T], nL2(M^)), which for every t > 
has unit mass and zero average. 

Denoting by h the right hand side of equation ()5.3|) we notice that /i > 0, for 
every g > 0, due to the truncation of the kernel. Indeed, 

(5.4) h = Qm,Mi9, g) + Mg> -g {g * \v\m,M) + Mg>0. 
Further, by analogy with Lemma HH] we can estimate Qm,M{g,g) as follows: 

(5.5) \\Qm,Mig,g)\\Li<CM\\g\\Li\\9\\Li, i<p<oo 

(there will be no loss of moments since the kernel Bm,M is bounded). Therefore, 

h e L°°([0,T],L^nL='(M^)), 

as soon as g is in the same space. The unique weak solution / G L°°{[0,T], Ll n 
L^(M^)) of ()5.3|) is then obtained from the following integral representation: 

(5.6) f{v, t) = fo{v) * E{v, t)+ [ h{v, r) * E{v, t - r) dr, 

Jo 

where * denotes the convolution in v, and E{v,t) is the fundamental solution of 

E(v t) = - e"^-*^* 

The regularity of / is then guaranteed by the classical parabolic regularity result 
Section 3.3], and we have the bound 



(5.7) ||/||i/2([0,T]xM^) < C'Af (||^||l2([0,T]xR'V) + ||/o||/fi(RiV))- 

We denote by T the operator that maps g into /. We next establish that for a 
certain choice of constants Ai and A2 this operator maps the set 



(5.^ 



B=\feL\[0,T]xR^)) f>0, [ fdv = l, [ fvdv = 0, 

Jrn J^n 

/ f\v\'^dv<Ai, / f dv < Al, for a. a. t e[0,T]] 



into itself. Indeed, the nonnegativity of / is evident from the integral representation 
()5.(jj) . since h > 0. The mass and momentum normalization conditions follow easily, 
since ioi g E B the collision term Qm,M{9, g) integrates to zero when multiplied by 
1 or V. It remains to verify the last two conditions in ()5.8p . 
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For the first of tliese conditions, multiplying the equation ()5.3p by \v\^ and inte- 
grating by parts we obtain: 



(5.9) 



^/ f\v\^dv + Ml f \v\'Uv < 2N + M I g\v\'Uv 



~^ / g g*\v\'^\v — v^:\m,M dv dv^ < 2N + {M — mk) / g\v\'^dv 
where k = €^{1 — 0^)/^. Therefore, taking g so that 

Jrn mk 



yields the differential inequality 

d 
di 

Then, by Gronwall's lemma 



[ f\v\^dv + M [ f\v\^dv<MA[. 

Jr^ Jri^ 



/ f Ivl"^ dv < max iA[, / /o|fp(if 
Jrn V Jjjjv 



Therefore, setting 



Ai = max iA[, / fo\v\'^dvj, 
^ Jrn ^ 

we obtain the required estimate. 

To obtain a bound of / in we integrate the equation against / and use the 
inequality to estimate Qm,M{g, g)'- 



It, I ^'^^+ / lv/prf^ + M / fdv 

<CM\\g\\L49\\L4f\\L^+M\\9\\L4f\\L^- 



(5.10) 

By Sobolev's embedding and interpolation, 

where < A < 1 is as in (j4.16j) . Therefore, dividing (j5.1(J|) by and taking into 

account that = ||/||li = 1, we get 

+ K4f\\%'~' + MII/IU. < {Cm + M)\\gh.. 



Using the inequality 



true for all e > 0, we find 



e 



^ll/IU^ + (1/5 + M)||/|U. <K, + {Cm + M)\\g\\L.. 
We then get by Gronwall's lemma: 

(5.11) \\f\\L2 <max{\\fo\\L^,P + ^\\g\\L^). 
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where 

(5-12) p = —, 7T and 7 = — — . 

^ ^ 1/e + M 1/e + M 

Choosing e < I/Cm we get 7 < 1. Therefore, we obtain the inequahty ||/||l2 < A2 
if we set 

A2 = max(||/o|U2,/?/(l -7)). 

It is straightforward to verify that the set B is convex and closed in the strong 
topology of Li([0,T] X M^), using Fatou's lemma and the fact that the second 
moment in |f | is uniformly bounded for g E B. Further, the uniform in time bounds 
assumed in the definition of B imply the continuity of Qm,M{g,g) in L^- We can 
then deduce easily that the solution operator T itself is continuous, based on the 
representation ()5.6|) . Finally, the bound for the second moment and the regularity 
estimate ()5.7p imply that the operator T maps B into its compact subset. By the 
Schauder theorem, this proves the existence of a fixed point for T in B, which is 
thereby a weak solution G L^{[0, T], n L^{R^)) of (Q. 

Our next goal is to pass to the limit as M ^ 00 and then as m ^ 0, to recover 
the solutions with the "hard sphere" collision kernel. To this end, we will show 
that the bounds set forth in the apriori estimates hold for the fixed point solutions, 
and are uniform in M (and m). First of all, using the computation ()5.9j) it is easy 
to conclude that the second moment is bounded uniformly in M, as soon as m > 0. 
Indeed, we obtain the following inequality for f = fmM, 



(5.13) ^/ f\v\^dv<2N-km[ 



f\v\^dv, 



so the required bound follows by Gronwall's lemma. 

Further, we see that for every m > 0, M > and for every T > 0, the solutions 
are in L°°{[0, T], Llp{R'^)), for every p > 1. To see this we take K > and introduce 
the truncated function 



xP, 0<x < K 

RP + pKP-\x - K), x>K. 



Then \l/p^i<-(x) is convex in x, continuously differentiable, and has a bounded second 
derivative. It also verifies conditions ()3.4j) - ()3.7|l . so Lemma 13.21 applies. Taking 
^p.ftrd'^n as a test function in the weak form of (|5.2p and arguing as in Lemma [3.41 



we get 

d 
di 



(5.14) 



/ /^p,ic(|t;|')+p(p-l)fcpm / / fMvl^Pdvdv, 

ilRJV J J{|i,|2 + |j,^|2<^} 

<2pAM I f\v\^dv I f \v\^P-^ dv + 2p{2p - 2 + N) [ f\v\^P-^dv. 
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Therefore, if we take 1 < p < 2, we can pass to the hmit as -ft" — oo in ()5.14j) using 
the monotonicity with respect to K and the bound of f m. L\. This imphes 

/eL-([0,T],LyM^)), 

for every 1 < p < 2, with bounds generally dependent on m and M. By induction, 
the same property is extended to every p > 0. 

We see further that the bounds in L\p are in fact independent on M. Indeed, 
estimating the middle term in ()5.14|) using the inequality 

\'V - V^\m,M < m + \v - 

and following the arguments of Lemma 13.41 we get 



(5.15) 



d 
It 



[ f \v\^P dv + p{p - l)kpm [ f\v\^Pdv<Ap[ f\v\dv I f \v\^'p dv 
+ (pm / f \v\^ dv + 2p{2p - 2 + N)) / f\v\'^P-^dv. 



This implies that for every T > fixed and every p > 0, the bounds of / = fm,M in 
L°°([0,T],L^p(M^)) are independent of M. 

Using the established Lip bounds and the fact that / G H^{[0,T] x M^) we can 
make rigorous the arguments of Lemma 14.71 and then proceed as in Lemmas 14.8114.91 
obtaining that 

feL^{[0,T],H^p{R'')), 

for every n = 1,2..., and every p > 0, with bounds independent on M. This will 
allow us to pass to the limit as M — > cxd in the weak form and to show that the limit 
solutions satisfy the equation with the kernel 

(m + \u\) b{u, a). 

We can then substitute the computation ()5.15|) by the argument of Lemma f3. 41 and 
find the bounds in L°°{[0,T], Llp{M.^)) that are independent on m and T. Arguing 
as above we can then pass to the limit as m — >^ 0. The limit solution obtained in 
this step will then satisfy the equation with the "hard sphere" kernel. 

Finally, in order to treat the problem with the initial data /o G L\f\L log L(M^) we 
can take a sequence /" G C^(M^) that converges to /o in L^(M^). Then, since the 
constants in the bounds for the energy and entropy from Section |21 are independent 
of ra, we can pass to the weak L^-limit in the equations. The fact that the bounds 
of the solutions are independent of T allows us to continue the obtained solutions 
to [0,2T], and by induction, to [0, oo). 

To study the regularity of solutions with initial data we use the parabolic 
regularity of the equation [2H] to find that / G i/H[0,T] x M^), for any T > 0. 
Using this fact in combination with the bound / G L°°([0, oo), we can 
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make rigorous the argument of Lemma 14.71 and then proceed as in Lemmas 14.81 and 
14.91 to find the infinite differentiabihty of the solutions. □ 

We now turn our attention to the steady equation ()().2|1 and give a proof of The- 
orem 15.21 One of the possible approaches consists of adapting the arguments devel- 
oped above for the time-dependent case. In fact, as a careful reader will easily check, 
practically all arguments in the above proof apply to the steady equation: the Gron- 
wall lemma arguments will be replaced by the inequalities obtained by dropping the 
time-derivative terms. The only point that would need more careful attention is the 
moment estimate ()5.15|) . which is not uniform in T. It can be replaced by a more 
elaborate argument for the moment bounds in the case of the truncated collision 
kernel. We will, however take another approach, which will allow us to obtain the 
existence of the steady problem as a consequence of the regularization properties of 
the time-dependent equation. 

Proof of Theorem \5.^A The proof of Theorem l5. ll enables us to construct a semigroup 
on the convex set C made of those functions in L\ fl L^(]R^) with unit mass and 
zero mean. Denote it by {St)t>o- Our bounds imply that for all t > 0, the range of 
St is compact in C. Therefore, for all n the equation 

fn = S2-nfn 

is solvable by Schauder's theorem. Since /„ = S'i/„, the sequence is contained in 
a fixed compact of C, namely Si{C). We can therefore extract a subsequence which 
converges towards some /. Now for all < n we have 

fn = 'S'2-fc/n 

(because 2"*^ is a multiple of 2~"'), and we can pass to the limit as ^ 00 using the 
continuity of the semigroup, thereby obtaining 

/ = for all > 0. 

Therefore f = Stf for all t which is a sum of inverse powers of 2. Since the set of 
such times forms a dense subset of M+ and since the semigroup is continuous with 
respect to t, we conclude that 

/ = Stf, for all t > 0. 

This ends the proof. □ 

6. Uniqueness by Gronwall's lemma 

We next show that under the assumption that the initial data has the moment 
of order 3 finite, the solution to the time- dependent problem is unique. The proof 
uses an argument based on a certain cancellation property of the collision operator 
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multiplied by sgn(/) [HIIII (see also [HI] for discussion). We show that this property 
yields the desired result for the operator with inelastic collisions as well. 

Theorem 6.1. Assume that /q G Ll{M.^); then the equation (jO.ip with the initial 
condition f{-,0) = fo has at most one solution. 

Proof. Assume that / and g are solutions of ()0.1|) . with the same initial data /q. Set 
h = f — g and H = f + g. Then h satisfies the equation 

(6.1) dth -A,h = ^ {Q{h, H) + Q{H, h)), 

with the homogeneous initial data. Now take a function ipe^x), a continuous ap- 
proximation of sgn{x). We can take 

f -1, x< -e 
ipeix) = < x/e, —e<x<e 
[ 1, X > e. 

Multiplying equation ()6.1|) by tpeih) (1 + |i;p) and integrating by parts we get 

^1 h^,{h){l + \v\^)dv + ^ I \Vh\\l + \v\^)dv 
M JfiN 2e J{\h\<e} 

-2N [ (P,{h) dv=l [ {Q{h, H) + Q{H, h)) ^,(/i) (1 + \v\'') dv, 



2 



where 



-X + X < —e 

(f>e{x) = I ipit) dt= { x"^ /2e, -e < X <e 

X — e/2, X > e. 

To estimate the right-hand side we can adapt the argument that is known to work 
in the case of the elastic Boltzmann equation (cf. PlEI)- Passing to the weak form 
we get: 

[ {Q{h,H) + Q{H,h))i;,{h){l + \v\^)dv 

^^•^^ = / / / H4hij.Sh'){i + \v'\') + hMK)i^ + K\') 

— hil)s{h) (1 + l^;!^) — h^sih^) (1 + |f*|^) } \u\ b{u, a) dadvdv^. 

Since + \v'J\^ < |f P + we can estimate the integrals of the first two terms 
in the braces as follows: 

{2 + \v\^ + \v^\^ — Cm — 7— 1^ — / / \h\H^{2+\v\^ + \v^\^)\v — v^\dvdv^. 

4 Jm J^N 

Subtracting the third term in the integral ()6.2|) and noticing that 

\h^P,ih)-\h\\<\h\Xeih), 



34 



I. M. GAMBA, V. PANFEROV AND C. VILLANI 



where Xe{x) is the characteristic function of the interval [—£,£:], we obtain the esti- 
mate for the first three terms: 



\h\ H^: {1 + \v^\ ) \v — v.^.\ dv dv^ 
+ / h{l + \v\'^)\v - v^\dvdv^. 

Jm J{\h\<e} 

The fourth term in ()6.2|) contributes with another integral like the first one above, 
so we finally get 

^ / hi/j,{h)il + \v\^)dv 
at J^N 

<2N (f)s{h)dv+ / \h\H,{l + \v^\'^)\v -v^\dvdv^ 
+^ / h{l + \v\'^)\v - v^\dvdv^,. 

2 Jrn J{\h\<e} 

Passing to the limit as e — > 0, we find: 

^ [ \h\ {l + \v\'^)dv 
dt J^N 

<2N [ \h\dv+ [ \h\{l + \v\^y/''dv [ H{l + \v\y^Uv 
Jrn Jrn J^n 

<C [ \h\ {l + \v\^)dv, 
Jrn 

since H is assumed to be bounded in LKW'^). Now since h{0,v) = it follows by 
Gronwall's lemma that h{t, v) = for all times. □ 

Remark. The uniqueness result of Theorem 16. II is most certainly suboptimal. We 
believe that the uniqueness could be obtained in the class of initial conditions with 
finite mass and energy, with no additional assumptions, similarly to the the classical 
Boltzmann equation [ST] . The main technical obstacle for such a result is extending 
the Povzner inequalities in the case inelastic collisions to the class of slowly growing 
piecewise linear functions ip studied in We believe that this can be overcome 
with a more careful analysis of the inelastic collision mechanism. 



7. Lower bounds with overpopulated high energy tails 

In this section we obtain pointwise lower estimates of solutions to ()0.1|) and ()0.2|) 
showing that the behavior of the high-energy tails of solutions is controlled below 
by "stretched Maxwellians" Aexp(— a|f p/^). The bounds are established by using 
the comparison principle based on the parabolic (elliptic) structure of the equations. 
The following proposition establishes the particular role played by the "stretched 
Maxwellians": they can be used as barrier functions in the comparison principle. 
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Proposition 7.1. Let g{v) be a nonnegative function with finite mass po = J^n 9 dv 
and moment of order one, pi = J^n 9 \v\dv. Then for every r > 0, and every K > 0, 
there is a constant a > such that the function 



,|3/2 



h{v) = ife""!''! 
satisfies 

(7.1) Ah-Q'{9,h)>0, for all \v\ > r, 
Further, choosing b > large enough, the function 

h{v,t) = Ke-''-<'+^^^'^"\ 

satisfies 

(7.2) -dth + Ah-Q-{9,h)>Q, 
for allt>0 and all v G M^. 

Proof. To prove inequality 1)7.11) we fix an r > 0, compute, 



9^2i„,i 3(2A^-1 

and use the estimate 



Ah=r-a'\v\- ''^''~'U \v\-'/']h 
' 4 I I 4 II 



Q {9, h) = h{g* \v\) < (pi + po\v\) h, 

to obtain 

(7.3) Ah - Q{g, /^) > ((^ - po)|^;| - Pi - ^^^^~^^ a\vr'/') h . 

If > Po, tlie factor on riglit-hand side of ()7.3p attains its minimum for |f| = r. 
Therefore, inequality ()7.H) holds for every a > a*, where a* is the positive root of 
the quadratic equation 

— a^- ^ a - (por + pi) = 

For the time-dependent operator, denoting by {v) = (1 + |f P)^, we obtain 

(7.4) Ah-dth= l^^\v\'{v)-\3a+{v)-'/^)-^{vr'/' + b^h. 

Choosing a so that | > po and then b > + Po + Pi we obtain inequality ()7.2|1 
and complete the proof. □ 



The established property of the function h{v) is used in next lemma to obtain a 
comparison result for the steady equation. 
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Lemma 7.2. Let f G Lj(R^) be a nonnegative smooth solution to fj0.2|) . with the 
mass po > 0. Then, there is a constant K > 0, such that 

(7.5) f{v) > Ke-^''\''\"\ 

for all V G M^, where a is a constant as in Proposition \7.1\ 

Proof. Assuming the smoothness of the solution to ()().2j) . there is a constant cq > 
and a ball B{vo, tq) with vq G and tq > 0, such that 

(7.6) f{v)>co>0, if veBivo,ro). 

The value of Cq (as well as tq and Vq) depend on the solution / and use the fact that 
Po > 0. 

Since equation ()().2|1 is translation invariant, we can take g{v) = f{v + vq)', then 

(7.7) Ag - Q-ig, g) = Ag - (g * \v\)g < 0. 

Applying Proposition l7.1l to the function g{v) with r = rg we find the barrier function 
h{v) = cq exp(— a|f I^/^), for which we have 

(7.8) Ah-Q-{g,h) = Ah-{g*\v\)h>0, for > tq. 
and 

g{v) > h{v), for |t>| < tq. 

Therefore, letting U{v) = g{v) — h{v), subtracting ()7.8p from ()7.7|) we obtain the 
inequality 

AU - {g*\v\)U <0, |t;|>ro. 

To prove that U{v) > everywhere we apply a form of a strong maximum principle 
(see, for example, P^) to the operator 

CU = AU -u{U + h)U. 

We can reduce the problem to proving that f/ > in a bounded domain. Indeed, 
the decay conditions on / imply that for every £ > we can find R > such that 
I [/(f) I < £ if |f I > -R. Then we have 

C{U + e) = CU -eu{g) <0, ro < \v\ < R 

and U + e > for \v\ = tq and |f | = R. The strong maximum principle then implies 
that U + e > for all tq < \v\ < R. Letting e go to zero we get 

U > 0, for all |f I > tq. 

In view of the inequality (|7.6|) this implies 

^?(^;)>Coe-'^H^'^ 

or, applied to the function f{v), 
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with K = coe~"'^°' . This completes the proof of the lemma. □ 

By using a version of the maximum principle for the parabolic operator, we obtain, 
in a similar fashion, the pointwise lower bound for the time dependent problem. 

Lemma 7.3. Let f G L°°([0, oo), Ll{M.^)) be a nonnegative smooth solution to fl().l|l 
with the initial data /q > Coexp(— aol^^P^^)- Then, there are positive constants K, b 
and a, generally depending on the solution, such that 

f{v,t)>Ke-''-^^''^'^\ 

for allt>0 and all v G M^. Further, if there is a constant Ci and a ball B{vq, tq), 
such that 

f{v,t) > ci, if V e 5(wo,ro), 
for all t, then the lower bound 

f{v,t)>Ke-^\'"\"\ 

holds uniformly in time, where now > 0, a > and b > will depend on ci, vq 
and tq. 

Proof. To prove the first statement of the lemma we use the second part of Propo- 
sition 17.11 and repeat the comparison arguments of Lemma 17.21 taking 

h{v, t) = if e-^*-"<">'^' 

and using the strong maximum principle for the parabolic operator on U = f{v + 

vo) - h{v) 

CU = AU- u{f) U - dtU. 

For the second part, the additional assumption made on / allows us to repeat the 
proof of Lemma 17.21 using the function h from ()7.4|1 . □ 

Remark. It is tempting to conjecture that solutions to ()().2|1 should satisfy a point- 
wise upper bound of the type i^'' exp(— a'|f |^/^), for certain values of a' and K'. 
However, the application of an argument based on the maximum principle requires 
estimating Q~^{f,f) pointwise, which is generally a difficult problem. Assuming a 
"no-cancellation" property in the spirit of the argument [87j . 

(7.9) Q+(/, /) - Q-(/, /) < -KQ'if, /), 

where fc^ > 0, a pointwise upper bound is indeed obtained by the maximum principle 
techniques. However, a justification of ()7.9|1 at the present time seems to be out 
of reach. Notice that quite recently Bobylev et al. [7j were able to prove an upper 
bound "in the sense" , namely, that for a certain choice of a' > 

/ f{v) exp(a'|t;|^/^) dv < +oo, 

which could possibly be a hint in favor of the pointwise bound hypothesis. 
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8. Concluding remarks 

We studied the existence, uniqueness and regularity for the time-dependent equa- 
tion (jn.lll and the existence, regularity for the steady equation ()0.2|) . An important 
problem that remained beyond the scope of our study is the convergence of the 
time-dependent solutions to the steady ones as time approaches infinity. In fact, 
this remains a serious open problem, since no Lyapunov functional for the time- 
evolution is known to exist. A number of other interesting questions can be raised 
in connection to the obtained results. Are the steady states unique up to a normal- 
ization? Do the steady solutions necessarily have radial symmetry? (This can be 
expected from the rotation invariance of the equations; the existence and regularity 
of radial solutions can be obtained by applying our analysis to the reduced one- 
dimensional problem, as in JT], or just by working in spaces of radially symmetric 
functions). 

We hope that the methods developed in the present work for the case of diffusion 
forcing could be useful for studying other problems involving the Boltzmann (En- 
skog) collision terms with other collision and driving mechanisms. In particular, a 
generalization to the case of a heat bath including a friction term seems to be rather 
straightforward. (The lower bounds in that case are expected to be Maxwellians.) 
It is also likely that applying the techniques of this paper should yield results for 
problems with the normal restitution coefficient dependent on the relative velocity 
I31IS1, which would allow us to study a broader range of physical phenomena. 

Another problem worth studying is the (quasi-)elastic limit a — > 1. The steady 
states for the Boltzmann equation with elastic interactions (a = 1) and vanishing 
diffusion (/i = 0) are Maxwellians, while for every fi > and every a < 1 we have 
a "3/2" lower bound. Obtaining quantitative information on the transition to the 
Maxwellian steady states would be valuable. We hope to address some of these 
questions in our future work. 
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